Experimental determination of the symmetry energy of a low density nuclear gas 
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Experimental analyses of moderate temperature nuclear gases produced in the violent collisions 
of 35 MeV/nucleon 64 Zn projectiles with 92 Mo and 197 Au target nuclei reveal a large degree of 
alpha particle clustering at low densities. For these gases, temperature and density dependent 
symmetry energy coefficients have been derived from isoscaling analyses of the yields of nuclei with 
A < 4. At densities of 0.01 to 0.05 times the ground state density of symmetric nuclear matter, the 
temperature and density dependent symmetry energies range from 9.03 to 13.6 MeV. This is much 
larger than those obtained in mean field Calculations and reflects the clusterization of low density 
nuclear matter. The results are in quite reasonable agreement with calculated values obtained with 
a recently proposed Virial Equation of State calculation. 

PACS numbers: 24.10.i,25.70.Gh 



In a recent theoretical paper, Horowitz and Schwenk 
have reported the development of a Virial Equation of 
State (VEOS) for low density nuclear matter This 
virial equation of state is a thermodynamically consis- 
tent equation of state, in which the virial coefficients for 
nucleon-nucleon, nucleon-alpha and alpha- alpha inter- 
actions are derived directly from experimental binding 
energies and scattering phase shifts. It is argued that 
contributions from 3 H and 3 He are expected to be small 
and they are ignored in the calculation. In the work re- 
ported in reference [l[ these virial coefficients were then 
used to make predictions for a variety of properties of 
nuclear matter over a range of density, temperature and 
composition. The authors view this virial equation of 
state, derived from experimental observables, as model- 
independent, and therefore a benchmark for all nuclear 
equations of state at low densities. Its importance in both 
nuclear physics and in the physics of the neutrino sphere 
in supernovae is discussed in the VEOS paper 0]. A par- 
ticularly important feature of the VEOS, emphasized in 
reference is the natural inclusion of clustering which 
leads to large symmetry energies at low baryon density. 

In this paper we extend our investigations of the nu- 
cleon and light cluster emission that occurs in near- Fermi 



energy heavy ion collisions @, 0, 0, H, HI to investigate 
the properties of the low density participant matter pro- 
duced in such collisions. The data provide experimental 
evidence for a large degree of alpha clustering in this 
low density matter, in agreement with theoretical pre- 
dictions [l[ 0, @, 0] ■ Temperature and density dependent 
symmetry free energies and symmetry energies have been 
determined at densities of 0.05po or less, where po is the 
ground state density of symmetric nuclear matter, by ap- 
plication of an isoscaling analysis I(J 11 1. The symme- 
try energy coefficient values obtained, 9.03 to 13.6 MeV, 
are much larger then those derived from effective interac- 
tions in mean field models. The values are in reasonable 
agreement with those calculated in the VEOS treatment 
of reference [l[ . 



EXPERIMENTAL PROCEDURES 

The reactions of 35A MeV 64 Zn projectiles with 92 Mo 
and 197 Au target nuclei were studied at the K-500 Super- 
Conducting Cyclotron at Texas A&M University, using 
the 4tt detector array NIMROD [jj. NIMROD consists 
of a 166 segment charged particle array set inside a ncu- 
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tron ball. The charged particle array is arranged in 12 
concentric rings around the beam axis. In those rings, 
the individual segments are fronted by ionization cham- 
bers (IC) filled with 30 Torr of CF 4 gas. Front and back 
windows were made of 2.0 /im aluminized Mylar foil. In 
each of these forward rings, two of the segments have 
two Si detectors( 150 and 500 /.tm thick) between the 
IC and Csl detectors (super telescopes) and three have 
one Si detector (300 thick). Each super telescope is 
further divided into two sections. The Csl detectors are 
10 cm thick Tl doped crystals read by photomultiplier 
tubes. For these detectors, a pulse shape discrimination 
method is employed to identify light particles. In all tele- 
scopes particles are identified in atomic number. In the 
super telescopes, all isotopes with atomic number Z=10 
are clearly identified. 



The energy calibration of the Si detectors was carried 
out using alpha particles from a 228 Th source and the 
observed punch through energies of identified reaction 
products. Since the energy losses of the lighter parti- 
cles, in particular the high energy Hydrogen isotopes, are 
rather small in the Si detectors, evaluation of the energy 
deposited in the Csl crystal from the energy loss in the 
Si detectors requires special care for higher energy par- 
ticles. Therefore, an additional energy calibration was 
performed employing Si detectors of thicknesses 1mm, 
backed by Csl detectors of three different lengths (1cm, 
3cm and 5cm) to measure the inclusive energy spectra of 
light charged particles from the reaction 64 Zn + 92 Mo at 
47A MeV The energy spectra were measured at all angles 
corresponding to those of the 12 rings of NIMROD. The 
combination of thicker Si E detectors and observation of 
high energy punch-through points for the particles which 
traversed these thinner Csl detectors allowed us to deter- 
mine the energy spectra with a high degree of confidence. 
We then used the 64 Zn + 92 Mo at 47A MeV as a stan- 
dard reaction to determine the Csl energy calibrations for 
all other runs. Neutron multiplicity was measured with 
the 4-7T neutron detector surrounding the charged parti- 
cle array. This detector, a neutron calorimeter filled with 
Gadolinium doped pseudocumene, consists of two hemi- 
spherical end caps and a cylindrical mid-section. The 
mid-section is divided into four separate 90 degree quad- 
rants. The hemispheres are 150 cm in diameter with 
beam pipe holes in the center and they are upstream and 
downstream of the charged particle array. Thermaliza- 
tion and capture of emitted neutrons in the ball leads to 
scintillation which is observed with phototubes providing 
event by event determinations of neutron multiplicity but 
little information on neutron energies and angular distri- 
butions. Further details on the detection system, energy 
calibrations and neutron ball efficiency may be found in 
reference Q. 



DATA ANALYSIS 

An inspection of the two dimensional arrays depict- 
ing the detected correlation between charged particle 
multiplicity and neutron multiplicity in NIMROD (not 
shown) , reveals a distinct correlation in which increasing 
charged particle multiplicity is associated with increas- 
ing neutron multiplicity. Although there are significant 
fluctuations reflecting both the competition between dif- 
ferent decay modes and the neutron detection efficiencies, 
these correlations provide a means for selecting the higher 
multiplicity, more violent collisions. For the analysis re- 
ported in this paper, we have selected events correspond- 
ing to the largest observed neutron and charged particle 
multiplicities. For 64 Zn+ 197 Au at 35 AMeV this selection 
corresponded to 442,000 events, while for 64 Zn+ 92 Mo at 
35 AMeV it corresponded to 280,000 events. Many of the 
techniques applied in this analysis have been discussed 
previously in greater detail in references 0, 0, [||. 
The detected correlation between charged particle mul- 
tiplicity and neutron multiplicity in NIMROD provides 
a means for selecting the collision violence 0, 0, H, 0]. 
For the systems studied, reference Q also contains an 
extensive comparison of earlier experimental results with 
predictions of Antisymmetrized Molecular Dynamics cal- 
culations (AMD-V) 

For the selected events we carried out analyses using 
three-source fits to the observed energy and angular dis- 
tributions of the light charged particles. The assumed 
sources are the projectile-like fragment source (PLF), the 
target-like fragment source (TLF) and an intermediate 
velocity (IV) source [H, [3, [H, [l6| . Even though the 
system evolves in a continuous fashion, such source fits 
provide a useful schematic picture of the emission pro- 
cess. From the fits we obtained parameters describing the 
ejectile spectra and multiplicities that can be associated 
to the three different sources. As in the earlier works, the 
IV source is found to have a source velocity very close to 
half of that of the projectile reflecting the initial decou- 
pling of the momentum sphere of the participant matter 
from that of the remaining nucleons. This important fea- 
ture of the dynamically evolving system manifests itself 
as kinematic differences between the early emitted light 
ejectiles (gas) and the remaining matter (liquid). In the 
following we probe the properties of this early ejectile gas 
as the system relaxes toward equilibrium and the two mo- 
mentum spheres become more and more similar. 

Treatment of the early cluster emission in a coales- 
cence framework [13, [l8| has proven to be a very effec- 
tive tool for understanding the energy spectra of the light 
ejectiles [13, LLJ]. We have previously developed exten- 
sions of these techniques to probe the early interaction 
zone and dynamic evolution in Fermi-energy collisions 
0, M, H, Si- For the present work we have derived 
information on the early thermalization stage of the re- 
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action 

H H M HI by focusing on the properties of 
early emitted mid-rapidity particles associated with the 
IV source. Such a selection minimizes contributions from 
the other sources. In addition, yields assigned to the 
TLF source are subtracted from the experimental yields. 
Thus, the yields of higher energy particles are relatively 
uncontaminated by later emission processes. AMD-V 
calculations reported previously [H, 0] indicate that the 
velocities of early emitted light particles decrease rapidly 
with increasing average emission time. We have exam- 
ined the evolution of the system by determining various 
parameters characterizing the ejectile yields, i.e., tem- 
perature, alpha mass fraction, N/Z ratios and isoscaling 
parameters, as a function of ejectile velocity. The ve- 
locity employed is the "surface velocity", V sur f, of the 
emitted particles in the IV frame, defined as the velocity 
of an emitted species at the nuclear surface, prior to ac- 
celeration in the Coulomb field [l3| . The energy prior to 
Coulomb acceleration is obtained in our analysis by sub- 
traction of the Coulomb barrier energy derived from the 
source fits. In earlier studies we have employed the cal- 
culated correlation from AMD-V predictions to calibrate 
the time-scale associated with our data 0, 0, [H, 0] ■ 

TEMPERATURE DETERMINATIONS 

To characterize the temperature evolution of the sys- 
tem we have used the yields of 2 H, 3 H, 3 He and 4 He clus- 
ters to determine the double isotope temperature, Thhc, 
as a function of V sur f in the IV frame [2l|, • 
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FIG. 1: Double isotope ratio temperatures, Thhb, as a func- 
tion of surface velocity for 35 MeV/nucleon 64 Zn + 92 Mo 
(open squares) and 64 Zn + 197 Au (open diamonds). Also 
shown is the average of the two (solid triangles) . Uncertain- 
ties are estimated to be < 30% 



to occur at times near 115 fm/c. Entry into the evapora- 
tion stage of the reaction occurs at V sur / near 3 cm/ns, 
corresponding to times near 160 fm/c. At that latter 
point, the temperatures are near 6 MeV, comparable to 
those observed in the caloric curve plateaus for these sys- 
tems [111 . In the plateau region average densities (liquid 
plus gas) as low as 0.4po are indicated by a previous ex- 
panding Fermi-gas analysis 24 1 . 
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where Ry sur /. , the double ratio of cluster yields, Y, for 
clusters with the same surface velocity, is [Y( 2 H) Y( 4 He)] 
/ [Y( 3 H) Y( 3 He)] and the constants 14.3 and 1.59 reflect 
binding energy, spin, masses and mass differences of the 
ejectiles. Equation [TJ for particles with the same V sur /, 

differs from the usual formulation by a factor of | 2 ap- 
pearing in the logarithmic term in the denominator Q . 

We present, in Figure [TJ results for the double iso- 
tope ratio temperatures as a function of surface velocity 
in the intermediate velocity frame. These temperatures 
have previously been reported in reference 4. The tem- 
perature evolution of the IV source ejectiles is seen to be 
quite similar for the two different systems. For each sys- 
tem investigated the double isotope ratio temperature de- 
termination exhibits a high maximum temperature, near 
13 MeV. In each case, after reaching a maximum the 
temperature then decreases monotonically as the veloc- 
ity decreases. In references we discuss this evolution 
and present varying pieces of evidence that at velocities 
smaller than the peak velocity, the system attains chem- 
ical equilibrium, at least on a local basis. This appears 



ALPHA MASS FRACTIONS 

The VEOS paper []J employs the alpha mass fraction, 
X^4, defined as the ratio of mass bound in alpha parti- 
cles to the total mass of the matter under consideration, 
to characterize the degree of alpha clustering at differ- 
ent densities and temperatures. We have determined ex- 
perimental values of as a function of velocity. For 
this purpose it was assumed that the unmeasured neu- 
tron multiplicity at a given velocity was the product of 
the 3 H/ 3 He yield ratio times the proton yield for that 
velocity. For the IV souce particles this is implicit in 
coalescence assumptions 03, EH, EH and consistent with 
experimental results 25j, [26| . In this way both the total 
mass yields and the mass fractions could be deduced. The 
sum of the mass emitted from the IV source is 31.9 ± 2.5 
amu for the Mo target and 29.7 ± 2.4 amu for the Au 
target. As in the case of the temperature evolution, the 
close correspondence of ejected mass provides evidence 
of production of a similar zone of participant matter in 
the two different cases. 

In Figure [2] we present the results for alpha mass frac- 
tions for the IV source ejectiles of both systems. For the 
two systems evolves in a similar fashion with surface 
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FIG. 2: Alpha mass fractions, Xa, as a function of surface 
velocity for 35 MeV/nucleon 64 Zn + 92 Mo (squares) and 64 Zn 
+ 197 Au (diamonds). Also shown by the dashed line is the 
average of the two. Solid triangles represent derived values 
used for comparison to VEOS results. See text. 



velocity. As the surface velocity decreases ~X.a increases 
dramatically. is smaller for the more neutron rich 
64 Zn + 197 Au entrance channel. Figures [1] and [2] clearly 
indicate that large alpha mass fractions are associated 
with the lower temperatures. 



ISOSCALING 

Horowitz and Schwenk have pointed out that extensive 
alpha clustering in the low density gas leads naturally 
to an increase in the symmetry energy for the cluster- 
ing system For comparison to the symmetry energy 
predictions of the VEOS model we have derived symme- 
try energies from the ejectile yield data by employing an 
isoscaling analysis. Such analyses have been reported in 
a number of recent papers [lCj, [ll 



, |2a 123, |30(. In 

this approach the yields of a particular species Y(N,Z) 
from two different equilibrated nuclear systems, 1 and 
2, of similar temperature but different neutron to pro- 
ton ratios, N/Z, are expected to be related through the 
isoscaling relationship 



Ce 



aN+pZ 



(3) 



where C is a constant and fJt(n) and fi(p) are the neu- 
tron and proton chemical potentials. 

The isoscaling parameters a — (^(n) — P L i{ n ))/T and 
= (a*2(p) — /xi(p))/T, representing the difference in 
chemical potential between the two systems, may be ex- 
tracted from suitable plots of yield ratios. Either param- 
eter may then be related to the symmetry free energy, 



Fsym- With the usual convention that system 2 is richer 
in neutrons than system 1, 



a = 4F s „ m ((Zi/Ai) 2 - (Z 2 /A 2 ) 2 )/T (4) 
= AFsy^Nx/Ax) 2 - (N 2 /A 2 f)/T (5) 

where Z is the atomic number and A is the mass num- 
ber of the emitter. Thus, F sym may be derived directly 
from determinations of system temperatures, Z/A ra- 
tios, and isoscaling parameters. We emphasize that the 
present analysis is carried out for light species character- 
istic of the nuclear gas rather than, as in most previous 
analyses, for the intermediate mass fragments thought to 
be characteristic of the nuclear liquid. 

In previous works attempts were made to employ sys- 
tems of similar Z and derive the a parameter which is 
expected to be less sensitive to Coulomb effects. It is 
normally assumed that all contributions to F sym other 
than that arising from the symmetry ene rgy e.g . bulk, 
Coulomb and surface, may be ignored 13, II, 2m 
r30j . Most commonly, the expression, 

a = 4 7 ((Z 1 /A 1 ) 2 - {Z 2 /A 2 f)/T 



28, 2 



(6) 



in which 7 is the symmetry energy coefficient, has been 
employed with estimates of Zj/Aj, Z2/A2 and T to ob- 
tain estimates of 7 [HI El EE 13 S Glfl]- While the 
approximate equation [5] may be reasonable for systems 
near normal density, it becomes increasingly less accurate 
at low densities where entropic effects become increas- 
ingly more important [3]]]. Indeed a recent paper has 
demonstrated the existence of isoscaling effects in a non- 
interacting system with no symmetry energy, resulting 
entirely from maximal entropy [3l[. Further, the VEOS 
treatment indicates that clustering in low density matter 
leads to entropic contributions to the free energy which 
can differ significantly in magnitude and, even in sign, 
from those obtained from mean field calculations In 
this work we employ equation [3] with experimentally de- 
termined isoscaling parameters, temperatures and Z/A 
ratios to determine the symmetry free energy, F sym . We 
believe the present work to be the first isoscaling analysis 
in which all relevant parameters of equations[4]and[5]have 
been determined experimentally. We then employ calcu- 
lated entropic corrections from reference [l[ to obtain the 
symmetry energy, E sym 



ISOSCALING PARAMETERS, a AND (3 

The isoscaling analysis of the two systems was done in 
the usual way by plotting isotope or isotone yield ratios, 
in our case for Z=l and Z=2 particles, and doing a global 
fit to the yield ratios to obtain a and j3. In Figure [3l 
the isoscaling parameters, a and 0, are plotted against 
surface velocity. 
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FIG. 3: Isoscaling parameters as a function of surface veloc- 
ity. The values for the parameters a (open circles), f3 (open 
squares and the parameter ratio a/ 13 (solid triangles). Esti- 
mated errors are ±5% 



FIG. 4: The quantity ((Zi/Ai) 2 - (Z 2 /A 2 ) 2 ) 2 as a function 
of surface velocity. Experimental determinations are repre- 
sented by solid circles. Results of AMD-V calculations are 
represented by open squares. 



We see that below 5.0 cm/ns the values vary slowly. 
Below 5 cm/ns the ratio, a//3, also varies slowly but in- 
creases dramatically at higher velocities, even changing 
sign. We have previously suggested that equilibration is 
observed for velocities below those for which the peak 
temperatures are seen. In our opinion the behavior of 
the individual parameters and of the a/ (3 ratio provides 
further evidence that equilibration is achieved for surface 
velocities below 5 cm/ns. 

Z/A RATIOS 

For each emitting system, the sums of Z and A for the 
emitted ejectiles, taken to be characteristic of that of the 
sampled gas, have been derived as a function of surface 
velocity. From these the values of ((Zi/Ai) 2 - (Z 2 /A 2 ) 2 ) 2 
were then obtained. These values are presented in Fig- 
ure [4] as a function of surface velocity. For comparison, 
theoretical values obtained for the gas fraction derived 
from the AMD-V calculation are also presented. The ex- 
perimental results and those for the gas in the dynamical 
AMD-V calculation are generally quite close in the 3 to 
7 cm/ns surface velocity. 

SYMMETRY FREE ENERGIES AND DENSITIES 

Experimentally determined values of T, Z/A and a 
were used in equation [4] to calculate F sym as a function 
of surface velocity. (Note that both temperature and 
density are varying with surface velocity.) The values of 
F sym are presented in column 3 of Table HI Values of 
range from 12.2 MeV to 5.4 MeV as the surface velocity 
decreases from 4.75 cm/ns to 0.75 cm/ns. 



While densities are not easily accessible experimental 
quantities in collision studies, knowledge of the densities 
at which the F sym determinations are being made is crit- 
ical to an interpretation of the measured values. 

As pointed out by Albergo et al, knowledge of the tem- 
perature allows the extraction of the free proton densities 
from the yield ratios of ejectiles which differ by one pro- 
ton, e.g., the yield ratio of 4 He to 3 He. Specifically , 

p p = 0.62 X l 36 T 3/2 e -19.8/T y ( 4jj e) jy (Sjj) (7) 

Here T is the temperature in MeV, Y refers to the 
yield of the species under consideration and p p has units 
of nucleons/cm 3 . 

Correspondingly, the free neutron densities may be ex- 
tracted from the yield ratios of ejectiles which differ by 
one neutron. Specifically, 

p n = 0.62 x 10 36 T 3 / 2 e - 20 - 6 / T r ( 4 He) /Y( 3 Re) (8) 

Once the free nucleon densities are known, the den- 
sities of the other particles may be calculated from the 
experimentally observed yields. This, again is done as 
a function of surface velocity. The results, obtained by 
summing the densities of particles with A = 1 to 4, for 
the two reaction systems, are presented in Figure [5] The 
values for the two systems show quite good agreement 
and indicate low densities. It is worth noting that our 
measurements of both the temperature and the associ- 
ated alpha mass fraction, X Q , also provide a means of 
estimating the densities by comparison with the Schwenk 
and Horowitz [l[ or Shen et al. § calculations. Like the 
Albergo calculation these calculations assume chemical 
equilibrium and lead to similar low densities. 
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FIG. 5: Estimated nuclear densities as a function of tempera- 
ture. Values derived from the experimental yield ratios using 
the Alberga formulas are indicated by solid squares. Theoreti- 
cal results from reference^ are shown by solid triangles. Val- 
ues derived from the theoretical calculations of reference |33| 
are shown for confining volumes of 5Vo (open squares) and 
10Vo (open diamonds). 



cles with A< 4, in reasonable agreement with that de- 
rived from the experiment. Although clusters are nat- 
urally produced in AMD-V calculations, early recogni- 
tion of clusters relies on a physical space coalescence ra- 
dius parameter. The value of 2.5 fm was used in this 
calculation. We then calculated, as a function of time, 
the average separation distance between all nucleons as- 
signed to the gas. From the average separation distance 
the average density is readily derived assuming spherical 
symmetry in the source frame. The AMD-V derived den- 
sities, averaged for Mo and Au targets, are in the range 
of2to4xl0~ 3 nucleons/fm 3 . The agreement between 
the AMD-V results and these experimental values in Fig- 
ure El while not perfect, is reasonable, particularly when 
it is recognized that the AMD-V results demand a sep- 
aration of nucleons for identification as "gas" while the 
clustering may actually occur in a slightly higher den- 
sity surface region of matter identified as "liquid" by the 
AMD-V analysis. Thus the results depend on the early 
recognition coalescence parameter. 



For the higher temperatures and densities, where the 
Shen et al. calculations include competition from heav- 
ier species, the Shen et al. calculation indicates slightly 
higher densities for a given temperature and alpha mass 
fraction. In the following analysis we employ the densities 
from the Albergo model treatment. The resultant den- 
sities, presented in column 4 of Table U and in Figure [5l 
range from 1.5 — 7.0 x 10~ 3 nucleons/fm 3 . For compari- 
son to these values we also show in Figure 03 theoretical 
results reported by Sil et al. [33| for hot Thomas-Fermi 
model calculations for 185 Re confined in volumes of 5 or 
10 times Vo, the nuclear volume at normal density. At 
temperatures of 5 to 10 MeV, they found average gas (or 
low density surface matter) densities in the range of 2 
to 10 x 10 -3 nucleons/fm 3 . Below T = 10 MeV these 
values are in good agreement with those derived from 
comparison of the data to the VEOS results and provide 
some support for those derived from the experiment. In 
an earlier application of the virial equation approach to 
similar types of reactions, Pratt et al. [34| derived densi- 
ties of 1 — 3 x 10 -3 nucleons/fm 3 and temperatures of 4.38 
MeV to 7.90 MeV. These are in the range of the present 
results. Additionally, it is worth noting that the derived 
gas densities at the temperatures in Table U are quite 
similar to those indicated by recent Fisher scaling model 
analyses of fragmentation data by Elliott et al. [351 ] . 

It is also possible to estimate gas densities using AMD- 
V model calculations once a calculational criterion capa- 
ble of distinguishing gas nucleons from liquid nucleons 
is established. For this calculation, we first verified that 
the parameterization of the AMD-V model reproduced, 
at 300 fm/c, a total mass fraction of gas, i.e. of parti- 



ENTROPIES AND SYMMETRY ENERGIES 

For the densities and temperatures derived from our 
data, we have employed the entropy equations of the 
VEOS calculation to determine the temperature and den- 
sity dependent entropic contribution to F sym . We have 
done this for the densities and temperatures of Table HI 

Given that the isoscaling parameter is sensitive to both 
the symmetry energy and symmetry entropy contribu- 
tions, derivation of the temperature and density depen- 
dent symmetry energy coefficient from the isoscaling pa- 
rameter requires an evaluation of the symmetry entropy 
contribution. We have based our extraction of this con- 
tribution on the entropy equations of the Schwenk and 
Horowitz paper. In that work, the entropy density is 
expressed as 



TABLE I: 



V surf 
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- 1 - sym 


P 




cm/ns 


MeV MeV 


nuc / fm 3 


MeV 


0.75 


3.65 


5.40 
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5P 

= — -n n log z n - n p log z p - n a log z a (9) 
+ TT + 2jX + ^ p z n b' pn ) 



T 



+ xs{z a b' a + 2z a (z n +z p )b' an ) 



(10) 
(11) 



Where P is the pressure, n„, n p and n a are the n, p and 
a-particle densities, z„, z p and alpha z Q are the particle 
fugacities, A(A Q ) denotes the nucleon (a particle) ther- 
mal wavelength and the bi terms denote the temperature 
derivatives of the second virial coefficients [l[ 

The entropy per nucleon is then 



S_ 
A 



s 

rib 



(12) 



where nf, is the total baryon density in the gas. Given 
that 3 H and 3 He are observed experimentally, we add 
their mixing entropy terms, contributions of the type, 
—riiliizi, to the entropy density equation. At the low 
apparent density of the nuclear gas with which we are 
dealing, we assume that the particle interaction terms 
are very small and that differences between them for the 
two different systems studied may be neglected. Thus 
we truncate the equation by dropping the interaction 
terms. We then replace the fugacities in the mixing en- 
tropy terms by fractional yields and calculate the entropy 
per nucleon for the two systems. With some scatter they 
exhibit an essentially linear dependence on surface veloc- 
ity. Thus, we have proceeded by making a linear fit to 
the entropy per nucleon as a function of surface velocity 
for each system. We then extract the entropy difference, 
AS, between the two systems. 

Knowing the total neutron and proton composition for 
the gas in the two systems we extract the symmetry en- 
tropy coefficient, S sym , 



S sym = AS/((N-Z)/A) 2 



(13) 



The values thus derived are presented as a function of 
surface velocity in Figure[6l For comparison we also show 
the values calculated from the Schwenk and Horowitz 
VEOS for the corresponding temperatures and densities. 
The agreement is quite reasonable except at the very low- 
est density. 

The symmetry energies obtained are presented in col- 
umn 5 of Table Q] The determination leads to values 
of the symmetry energy coefficients which range from 
9.03 to 13.6 MeV. We estimate the overall errors of these 
to be ±15%. The symmetry energy coefficients are also 
plotted against density in Figure [7] where they are com- 
pared to those of uniform nuclear matter which are pre- 
dicted by the Gogny effective interaction [36j and to the 
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FIG. 6: Ssym vs v sur f. Experimental data is represented 
by closed squares and calculated values are shown by open 
circles. 




FIG. 7: Derived symmetry energy coefficients as a function of 
baryon density. Solid diamonds indicate results using densi- 
ties of column 4 in Table [T] Solid line indicates the variation 
predicted by the Gogny interaction. The dotted line repre- 
sents the function 31.6(p/p ) 105 



31.6(p//9o) 105 dependence suggested by a recent analy- 
sis of isospin diffusion data [371 ] . The derived values of 
E S j, m are much higher than those predicted by mean field 
calculations which ignore the cluster formation. 



CONCLUSIONS 

For nuclear gases with average proton fraction, Y p ~ 
0.44, corresponding to the experimental average asym- 
metry of the Mo and Au target results and densities at 
and below 0.05 times normal nuclear density, experimen- 
tal symmetry energy coefficients of 9.03 to 13.6 MeV, 
with an estimated uncertainty of 15% have been derived 
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from experimentally determined symmetry free energies, 
F 8ym , determined using the isoscaling method. The sym- 
metry energies are far above those obtained in common 
effective interaction calculations and reflect cluster for- 
mation, primarily of alpha particles, not included in such 
calculations. The entropic contributions to the symme- 
try free energies derived from isoscaling analyses are very 
important at low densities and are strongly affected by 
the cluster formation. The experimental values obtained 
agree well with the predictions of the Virial Equation 
of State model proposed by Horowitz and Schwenk 
Inclusion of 3 H, 3 He and heavier clusters into the theo- 
retical formalism is certainly desirable. Measurements at 
slightly higher densities should be more sensitive to the 
interaction terms in the VEOS. 
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